

c= 21839262149008049274877780752446238293275890311566388523279198938849808962703425641238786067467508703093919234523232663053109522714732976531162075009707830309860135771109239220561092640496155709721227114461407454500057862947345932212052683967460416324819852661346135058013412441373992690839146716101809134598
n= 68169144741149098870815008157932118480267780082160027894190585231761354416244478110260464767604477038264129293569208701810025783427850185954992631235591973586615738638203024593579825895524769875068215844631524125282275486083378947501116174414797229145722176836098299961083840770249199297282599562339498107737
e=65537

import libnum
import gmpy2

p1,s=gmpy2.iroot(n,2)
q=gmpy2.next_prime(p1)
p=n//q
print(n==p*q)
phi=(p-1)*(q-1)
d=libnum.invmod(e,phi)
m1=pow(c,d,n)
#m1=m*(p-2)!%p
#(p-2)!%p=1
#m1=m
print(m1)
print(libnum.n2s(int(m1)))
